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Abstract 

Using a simple cellular automaton with stochastic rules we show 
the possible emergence of thermally activated avalanches (power law 
distributed) in type-II superconductors. Scaling relations between 
the exponents characterizing these distributions and those obtained 
from field driven experiments are derived and proved through simula- 
tions. It is also shown that the conditions for the appearance of these 
avalanches are independent of the pinning mechanism. The relevance 
of our simulations for recently reported experimental results is also 
outlined. 

PACS: 64.60.Lx, 74.69.Ge 
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1 Introduction 

Magnetic flux penetrates type-II superconductors above a certain critical 
field H c i in the form of vortices. The interaction of these vortices with the 
pinning centers produces a magnetic flux profile inside the superconductor 
with a slope proportional to the critical current density inside the sample, 
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j c , defined as the maximum current density the material supports without 
dissipation, a situation which is accounted by the so called Bean's Critical 
State Model [[!]. This picture, as de Gennes fristly noted [0], is very similar 
to the case of sandpiles, where a constant slope appears in the pile resulting 
from the competition between gravity and the friction between grains. 

In 1987, Bak et al |3|, [|] proposed a theory - now known as Self Organized 
Criticallity theory (SOC) - to explain the existence of self similar structures 
in Nature. Since then, SOC has been used to interpret the dynamics of 
many size avalanches in sandpiles^, earthquakes [§], evolution |7|], and other 
phenomena, see || for a general review. 

The ocurrence of self-organized criticallity was soon searched also in su- 
perconductors where field driven experiments have been designed [[], [IT], [TIJ 



and many numerical simulations developed [lj| [L3|, [14], [15], [16], |17|, unfortu- 
nately without conclusive answers. 

Superconductors differ from most systems exhibiting SOC by the rele- 
vant role played by the temperature. The temperature causes the relaxation 
of the critical state leading to a nearly logarithmic magnetization decay, 
m(t) ~ ln[t) ||18|| . In the early 90's many researchers tried to relate the 



role played by the temperature to the existence of many size avalanches in 



relaxation experiments [19], |2CJ, |21], p23 |. However in 1995 Bonabeau and 



Lederer [ 21 , 25|, approximately solved the difusion equation for the mag- 
netic field inside a superconducting slab and demostrated that (within the 
usually accesible time scales in the experiments) it is impossible to determine 
the existence of thermally activated avalanches by classical magnetic relax- 
ation measurements, i.e by the study of the decay of the mean value of the 
magnetization in the sample flT8|, |23|| . 

In 1998 Aegerter studied the magnetic relaxation of a Bismuth sin- 
gle crystal but instead of the usually measured mean value of m(t) |TB[ he 
put attention to the fluctuations during the decay of the magnetization and 
showed evidences of power law distributed thermally activated avalanches. 

In this work we develop a simple scenario able to account for the existence 
of these thermally activated avalanches. This does not mean we claim for 
the existence of SOC during the relaxation of the magnetization. SOC is 
well defined only for a system in a marginal stationary state and it is not the 
case for the vortex lattice in the presence of thermal activation. What we are 
claimining is that, because of the complex interaction between vortices, the 
pinning centers and the temperature, many size avalanches of moving vortices 
may produce the relaxation of the critical state as previously determined in 



2 



reference IOB 



The remaining of the paper is organized as follows. In the next section 
we describe the Cellular Automaton used in our simulations. In section 3 
we present and discuss the numerical results. Then, section 4 is devoted to 
the study of the scaling relations between our distributions and those usually 
obtained in field driven experiments. In section 5 we outlined some conditions 
needed for the ocurrence of many sizes thermally activated avalanches and 
finally in section 6 the conclusions are given. 

2 The model 

While the use of "real" forces between vortices in molecular dynamics simula- 
tions 1T3, O better resemble the experimental situation than simple cellular 



automata, they are by far more time consumming and it is an important 
drawback of the method, specially when we are looking for critical exponents 
or when we introduce the effects of the temperature in the system. 

Recently, to neglects these problems Bassler and Paczuski (jTJJ] introduced 



a simple cellular automaton to study the behaviour of the vortex lattice in 
type-II superconductors. This cellular automaton avoids part of the relevant 
physics of the vortex lattice such as the variation of the pinning strenght 
with the increasing field, the possible mistmach between the vortex lattice 
and the pinning centers, the elasticity of the vortex lattice, etc. However, 
it contains the interaction between vortices and with the pinning centers, 
the long range order of the vortex interaction, first by the introduction of 
the parameter r (see below), but also implicity assuming that each lattice 
cell contains more than one vortex. In addition it is able to predict the 
self organization of the lattice in a critical state characterized by power law 



distributed avalanches ||14|, |15| and the irreversibility of the magnetization. 
Then, aimed at describing the influence of the temperature on this critical 
state we adopt this model with some modifications. 

The cellular automata consists in a 2D honeycomb lattice, where each 
site is characterized by the number of vortices on it m(x), and by its pinning 
strength V(x) equal to with probability p, and to q with probability 1 — p. 
The force acting on a vortex at site x in the direction to site y is calculated 
as: 

F x ^ y = -V(x) + V(y) + (m(x) - m{y) - 1) 



3 



+r[m(x 1 ) + m(x 2 ) — m(y 1 ) - m(y 2 )] (1) 



where X\ and x 2 are the nearest neighbors of x (other than y) while yi and 
y 2 are the nearest neighbors of y (other than x) and r is a measure of their 
contribution to the total force on the vortex x (0 < r < 1). A vortex in the 
site x, moves to its neighbor site y if the force acting on it in that direction 
is greater than zero. If the force in more than one direction is greater than 



zero, then one of them is chosen at random [14] , |15[ jig] . 

To introduce the effect of the temperature we assumed that sites where 
the forces are lower than zero still have a probability of motion given by: 

P x ^ y ~ exp(-f/(j) AT) (2) 

where k is the Boltzman's constant and T is the temperature. The current, j, 
was locally calculated using the gradient of m(x) and U(j) represents different 
pinning barriers proposed in the literature U(j) = U a j c /j, U(j) = U \n(j c /j) 

andU{j) = U {l-j/j c )m- 

An avalanche starts by randomly choosing a lattice site, and calculating 

(fj). If it is smaller than a random number the procedure is repeated, else the 

vortex moves perturbing its neighbours. Then, the direction of motion of the 

new unstable vortices is calculated using (p. At this point, all the sites are 

updated in parallel until no more unstable sites persist. The avalanche size is 

defined as the number of topplings corresponding to the thermal activation of 

one vortex while the avalanche duration is defined as the number of updatings 

necessary to complete one avalanche. 

In all the cases the procedure was repeated for 10 4 m.c.s, were one m.c.s 
was defined by the L 2 calculation of (|2|), and lattices up to L = 200 were 
used. The initial configuration was obtained slowly adding vortices to the 
system (at T = 0) until a critical slope was reached [pf] . The boundaries 
"parallel to the net vortex motion" were assumed periodic while the other 
two were fixed to mimic the applied external field. 

The magnetization, M was calculated as the mean magnetic field inside 
the sample minus the external applied field |T7|], i.e., 



M = l ^B(z)-H (3) 

i=0 

where H is the field, i.e. the number of vortices, at the borders of the lattice. 
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3 Numerical results 



Figure 1 shows typical relaxation curves obtained for systems of different 
sizes using a vortex glass- like potential U(j) ~ j c /j |27fl and the algorithm 
described above. In Figure 2 is represented the relaxation curve for a sys- 
tem with T ~ oo, it means when we disregard the avalanche- like behavior 
previosly explained avoiding the calculation of equation ([!]) after a thermal 
jump. 

In both figures (see also the inset) three regimes are present, a plateau, 
then a logarithmic relaxation, and finally another plateau due to finite size 
effects. Only the time scales for these regimes are different, but this is irrele- 
vant from the experimental point of view. So, as already noted before |24|. |25 



our results suggest that it is not possible to decide about the existence or 
not of thermally activated vortex avalanches from "simple thermodynamic 
magnetic" relaxation measurements. Other pinning potential as well as dif- 
ferent Uo/kT relations were used |27| and no fundamental differences with 



the previous results were obtained. 

Figure 3 and 4 represent the integrated^ distribution of avalanche sizes, 
Di n t(s), and the integrated distribution of avalanche times, Di nt (t) obtained 
using a classical Anderson-Kim potential, U(j) = U {1 — j/j c ) |]2SL ETfl, for 
a system with L = 200 and U Q /kT = 10, as before, other pinning potentials 
were also used, resulting in a similar behavior. These distributions were 
obtained using the avalanche sizes and times (defined in section 2) obtained 
during all the relaxation process. 

As figure 3 an 4 clearly show, many size avalanches emerge. It does not 
mean the system is critical, instead it is relaxing from a critical state to 
its corresponding thermodynamic equilibrium. What we are showing is that 
this relaxation could proceed by means of many size avalanches in accordance 
with recent experimental results However, somehow more surprisingly, 
we will show in the next section that the exponents characterizing these 
distributions are related through simple scaling relations to the exponents 



derived in the context of SOC for systems in a critical state JR], [15 
Considering that D int (s) follows a power law: 

D mt {s) ~ s~ T - (4) 
the estimated exponent form figure 3 was r n = 2.70 ± 0.1 (different form the 
1 The meaning of this name will be clarified below 
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t = 1.63 obtained in references [I4"| and [Hi] for a field driven experiment) 
and, assuming D int {t) ~ t~ rtn for the integrated distribution of avalanche 
times, we got from figure 4, r tn = 4.0 ± 0.2. 

It is worth to mention here that the exponent r n was also experimentally 
determined in reference |26| and reported as r n = 2.0, lower than our value. 
This divergence can be explained since figures 3 and 4 represent the distribu- 
tion of avalanches obtained for all the relaxation process, i.e starting at the 
critical state and finishing at equilibrium, a situation imposible to account 
for in real experimental situations. 

To determine the distribution of avalanche sizes, using just part of the 
relaxation curves, gives different numerical estimates for r n and Tf n as is 
evident from figure 5. In fact, figure 5 represents five avalanche size distri- 
butions, P(s), obtained for different time intervals of the relaxation curve, 
from the upper to the lower curve, t — 1 — 10, t — 11 — 100, t = 101 — 1000, 
t = 1001 — 10000 and t = 10001 — 100000 m.c.s, which superposition corre- 
sponds to the full relaxation of the system (see figure 1). The straight line 
represents the integrated distribution of avalanche sizes, D int (t), obtained in 
figure 3, r = 2.7. Then, from the figure we can conclude that can be pre- 
dicted different exponents depending on the range of time measured. For low 
enough times, the exponent is lower than that associate to Di nt , while for 
large times a peaked distribution is obtained with only very small avalanches. 

Another source for discrepancies between our numerical estimates and 
experimental situations comes from the change of regimes of relaxation. In 
fact, there is not a priori justification to assume that many size avalanches will 
dominate the relaxation process within the all range of j and T, a situation 
that was deeply analyze in references |24| , |25| and is discussed in a different 
context in section 5. 



4 Scaling relations 

Rather than to introduce directly the derivation of our scaling relations we 
prefer to start with a short review of some important scaling concepts of the 
theory of self organized criticallity. 

Following the first ideas of Bak and collaborators 0, |], those systems 
who behave as predicted by SOC show a distribution of avalanche sizes and 
times that follow power laws, i.e , 
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and 



P(s) ~ s- T 



(5) 



P(t) ~ t~ rt (6) 

respectively. For systems not exactly in the critical state, these expresions 
transform into: 

P(s) ~ s- T /(s/5c) (7) 

and 

P(t)~r*f(t/t e ) (8) 

where s c and £ c reflect the departure of the system from criticality, s c ~ 
(jc — j) 1 ^ CT1 and ~ (jc — j) _ 2 being <7i and o"2 new critical exponents 
and where the function f(x) has the following properties f(x) — > cte if a: — ^ 
and f(x) —>■ if x — > cxd in order to recover the "critical picture" when j ~ j c . 

In finite size systems s c and t c also reflect the effect of the sample dimen- 
sions through two new critical exponents D and z. In fact, in analogy with 
the theory of critical phenomena s c ~ L D and t c ~ L z . Furthermore, the 
coherence length £ diverges at the critical state as: 

Z~(jc-j)-" (9) 

From the definitions of s c , t c and the equation (^) it is straighforward to 
show that s c ~ £ 1/uai and t c ~ ^ l / utJ \ Moreover, since for a finite size system 
at the critical state £ = L, our first scaling relation takes the form: 

- = ^ (10) 

Z cr 2 

As already discussed above, the integrated distributions of avalanche sizes 
and times, calculated in section 2, D int (s) and D int {t) result, since the system 
is relaxing, from avalanches obtained for values of current densities ranging 
from j c to j. These distributions are different from those obtained in typical 
field driven experiments or simulations since the last are obtained "in prin- 
ciple" for a fixed value of current density j c which, indeed, determines the 
criticality of the system. 
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Then, it is natural to assume that D int (s) and D int (t) are related to 
the distributions obtained just at the critical state, D(s) and D(t), by the 
following formulae: 



An*(s)~ I" S- T f(s/S c )dj (11) 
J 3c 



and 



Ant(*)~ /° S- Tt f{t/t c )dj (12) 

which inmediately explain the meaning of the label "integrated" used for 
these distributions. 

Then substituing the definitions of s c and t c in equations (0) and (p~2|) 
and after a simple change of variables, we obtain the following expresions for 
the integrated distributions of avalanche sizes and times: 

Ant(s) = s ~ T+ai f ( ' 1 crix ai ~ 1 f(x)dx (13) 



o 



D int (t) = s- Tt+ ^ / a 2 x^- l f(x)dx (14) 



o 

which prove that for s large enough both integrals are constants, and there 
is not a cut-off length in the integrated distributions, result already obtained 
in our simulations (see figures 3 and 4). 

Also from equations ( |T3"D and flT4] ) and the definitions of D int (s) and 
Di n t(t) we can inmediatly obtain the following scaling relations 

T n = T + (Ti (15) 



Ttn = r t + 02 

which in combination with ( |ToD leads to: 

-D Tti — t 



Tin - Tt 



(16) 



(17) 



In this way expresion (|T7|) establishes a connection between the exponents 
obtained in field driven experiments or simulations, T,r t ,D,z and those from 
thermally activated avalanches T n ,r tn . In fact, the results obtained in our 
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simulations and those obtained in references 1*3, IS] hold the previous rela- 
tion. 

However, some points deserve further discussion. The power law diver- 
gence of s c , t c and £ are strictly valid close to the critical state, j c . Far away 
from this state these divergencies no longer exactly hold, however considering 
the good results obtained in the check of our calculations and the scaling law 
(|T7| ) we believe that this last assumption is not relevant for the solution of 
the model. Also, the scaling law (|17|) was obtained assuming the complete 
relaxation of the system, so it is difficult to be proved in real experiments. 



5 Applicability 

Our previous picture assumes that a thermally activated vortex jump would 
affect its neighborhood generating an unstability that leads to a cascade of 
vortex jumps related to the vortex distribution into the sample. 

However, it is well known the existence of a characteristic time for ther- 
mally activated phenomena = t D exp (U(j)/kT) representing the time a 



vortex spends at a pinning site before jumping due to thermal activation |2? 

This means that our model will be valid if these avalanches occur within 
times lower than t t h, i-e. the avalanches should develop fast enough to be 
mutually independent. This resembles the idea developed by Vespignani et 
al in the context of sandpile and forest fire models. They showed through 
simulations and mean field considerations that one necessary condition for 
the occurrence of SOC, at least in these models, is the separation of time 
scales between the external exitation and the response of the system. 

Then, as mentioned above, the maximum time an avalanche persits is, 
tc — ^co(l — j / jc)^ 1 ^ 2 where t co is the time a vortex spends moving from one 
site to another, and of course depends on the local current and flux density 
in the system. Considering that the vortices are separated a distance a the 
time they spend traveling this distance is 

a 

tea =~ 18 

V 

where v depends on the Lorentz force acting on the vortex v = j^? /r], and 
a = (Qo/B) 1 / 2 which inmediatly gives the following dependence of t co with j 
and B. 
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In the critical state B varies along the sample. This variation is, even in 
the presence of thermal activation, very well accounted by the Bean model |3], 
p0[ . This means that for a fully penetrated sample 

B(x) = fM H - Hojx (20) 

where H is the external field. Then, substituing equations (|19D and fl20D in 
the definition of t c we obtain that 

*.= (21) 
j$o' VHoH - /j, jx 

Now, to determine the regime of applicability of our model, we must verify 
under which conditions the inequality t c « t t h holds. 

From the experimental point of view, the relevant avalanches to be de- 
tected measuring the fluctuations in the magnetization decay [^6| are those 



starting at the border of the sample, since are those who produce changes in 
m. Moreover, the avalanches starting at the border of the sample are also 
those with larger duration times since they have a larger area for spreading, 
(remember that the critical state in a type-II superconductor is symmetric 



with respect to the center of the sample). Then, we may assume in (21) 
x = and obtain the following inequality: 

^yl 3 i^ 2 « t exp(U(j)/kT) (22) 
which can be written as: 

J (l-j/j c )iM ex ^)l kT ) ( 23 ) 

where j* = ^f§t Q . Assuming, for example U(j) = U \n(j c /j) |2jJ, the previ- 
ous inequality takes the form: 

3 -^o--m~ i/n «i (24) 

Jc 

a = U /kT. 
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It is then straighforward to demostrate that equation fl2"^D holds under 
the following conditions, if a » 1 j must be much lower than j c (j << j c ). 
In the opposite case j* << j << j c . Similar expresions can be derived for 
the Anderson-Kim potentials and from potentials derived from the Collective 
Pinning Theory |27f . 

These conditions are the consequence of the competition between the 
increase of tth when j — > and the divergence of t c when j — > and j — > j c , 
see equation (|21|), and can be interpreted in the following way. Close to j c the 
avalanche durations are very high because the avalanche sizes become huge, 
so one always need to be far from j c , a situation often accounted in high 
temperature superconductors ljT8| , to assure that t av « t th . Particularly for 
U « kT, when thermally activated jumps become frequent (t t h small), high 
enough currents (j » j*) are also neccesary to assure rapid vortex motion 
during the avalanche, and hence low avalanche time durations. From the 
experimental point of view these conditions should be seen with some caution. 



For example, since j c decays with temperature ||27|| , for U a « kT, the range 
of current densities where thermally activated avalanches could appear is 
still narrower than that suggested by a simple inspection to the formula 
j* << j << j c so, we strongly recommend to look for these avalanches at 
low temperatures and in very disordered systems were U Q » kT. 

In the light of these results it is useful to come back to the experiment of 



Aegerter [pqj . He found one critical exponent characterizing the avalanche 
size distribution during the relaxation of the magnetization and that this 
exponent was independent of the temperature of the system. Neither of 
these results contradict our model. Even when his critical exponent was 2.0 
and our r n = 2.7, figure 5 indicates that small exponents are associated to 
small relaxation times in our model. This suggest that, if in the experiment of 
reference PEJ the time window had been shifted to larger times, and exponent 
closer to our one would have been observed. This does not mean, of course, 
that such a shift can be trivially performed in the practice. 

In addition he found, during the relaxation, one initial regime where 
avalanches are not power law distributed. While he explained this due to 
a transient period the system takes to reach the SOC, our results suggest a 
different explanation. During this period the system is still too close to the 
critical state j ~ j c , and power law avalanches are not yet developed since 
the thermally activated avalanches overlap each other. This explanation is 
consistent with the long time associated to this transient period, and to 
the dependence of this time with the temperature. Experimentally Aegerter 
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found that larger times are associated with larger temperatures and in fact, 
in our model larger temperatures imply the neccesity of lower values of the 
relation j/j c to found power law distributed avalanches and this means larger 
transient periods. 

6 Conclusions 

In conclusion, we developed a simple scenario to explain recently reported 
thermally activated avalanches power law distributed for type-II supercon- 
ductors. We proved that the exponents associated with these distributions 
depend on the time interval of the measurement. We also proved that the 
exponents characterizing a distribution of thermally activated avalanches ob- 
tained during the whole relaxation experiment, (i.e, from the critical to the 
equilibrium states), are related to those obtained in field driven experiments 
by scaling relations, a situation also supported by our simulations. The con- 
ditions for the appearance of these avalanches were discussed and it was also 
proved that, in a rough approximation, they do not depend on the pinning 
mechanism in the sample. All ours theoretical predictions are consistent with 
the known experimental results. 
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Figure captions 



Figure 1 Magnetic relaxation curves for systems of sizes L = 60, 100, 200. 
U(j) ~ jc/j, U /kT = oo. The inset shows the data collapse of the curves. 

Figure 2 Magnetic relaxation curves for systems of sizes L = 60, 100, 200. 
U(j) ~ jc/j, U /kT = 10 The inset shows the collapse of the curves. 

Figure 3 Avalanche sizes distribution for L = 200, U(J) ~ (1 — j/j c ), 
U /kT = 10. 

Figure 4 Avalanche times distribution for L = 200, U(j) ~ (1 — j/j c ), 
U /kT = 10. 

Figure 5 Avalanche size distribution for L = 200, U(j) ~ jc/j-, U Q /kT = 
10. From the upper to the lower curve: t = 1-10, t = 11-100, t = 101-1000, 
t = 1001 - 10000 and t = 10001 - 100000 m.c.s. The straigh line represents 
a power law with exponent 2.7. 
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